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Abstract 

We present models and computational results which indicate that the spatial and 
temporal regularity seen in Proteus mirabilis swarm-colony development is largely an 
expression of a sharp age of dedifferentiation in the cell cycle from motile swarmer 
cells to immotile dividing cells (also called swimmer or vegetative cells.) This contrasts 
strongly with reaction-diffusion models of Proteus behavior that ignore or average out 
the age structure of the cell population and instead use only density-dependent mech- 
anisms. We argue the necessity of retaining the explicit age structure, and suggest 
experiments that may help determine the underlying mechanisms empirically. Con- 
sequently, we advocate Proteus as a model organism for a multiscale understanding 
of how and to what extent the life cycle of individual cells affects the macroscopic 
behavior of a biological system. 

Key words: Proteus mirabilis, sweirm colony, age structure, space structure, par- 
tial differential equations, cell cycle . 

1 Introduction 

We present structured multiscale models and computational results of Proteus mirabilis 
swarm-colony development which indicate that the spatial and temporal regularity of the 
colony development is largely due to the regularity of the cell cycle, namely a sharp age of 
dedifferentiation from motile swarmer cells with a multinuclear structure to single-nucleus 
dividing cells (also referred to as swimmer cells or vegetative cells elsewhere in the literature) . 
We refer to our models as "structured multiscale" because they link the cellular scale to the 
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colony scale, and represent the age dynamics of individual cells using a population structure 
variable for age and an explicit aging term in the differential equation. 

This explanatory mechanism precludes the use of reaction-diffusion models that either 
average over the age variable, such as the models by Medvedev, Kaper, and Kopell [21], or 
ignore the age structure from the outset, such as the models by Czirok, Matsushita, and 
Vicsek m. 

The necessity of explicitly retaining the cell cycle through an age variable, rather than 
using density-dependent processes alone, makes Proteus an ideal organism for a multiscale 
understanding of how the kinetics at the individual cell level determines the macroscopic 
behavior of a biological system. 

The models in this paper are based upon those presented in [^j. The previous models 
assumed a sharp age of dedifferentiation between motile swarmer cells and immotile dividing 
cells. This was incorporated into the model using an upper bound on the age domain and 
then assuming that the swarmer-cell population density was zero above this maximum age. 

The models in j2] used the critical insight of Esipov and Shapiro in JU] that the age 
dynamics of individual swarmer cells underly the observed colony regularity and should be 
incorporated into the model explicitly, using a population structure variable in the partial 
differential equation for the swarmer-cell population. The models in differed substantially 
from those in ITUj by using a simpler diffusivity which depended only on the swarmer-cell 
density rather than an elaborate and unnecessary memory field - one that proved to be 
insufficient in generating the desired colony behavior when accurate numerical methods were 
applied to the model equations. The models in [2] also incorporated the known mechanism 
of a density-dependent lag in differentiation between dividing cells and swarmer cells. Other 
differences included the use of a radially-symmetric geometry rather than a linear geometry, 
diffusion based on isotropic random motion rather than Fickian (symmetric) diffusion, and 
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accurate numerical methods. The models in (TUI, when solved accurately, can generate a 
regular temporal cycle and contain a control of the ratio of swarm time to consolidation time, 
but are unable to generate the regularity of the spatial structure, namely the bulls-eye pattern 
with equally spaced concentric terraces. An appendix in j2] contains a detailed discussion 
and computations of the Esipov-Shapiro model, and the manuscript as a whole address in 
more detail than here the differences between the two modeling approaches. The focus in 
this paper is the need to represent the age structure of the swarmer-cell population explicitly, 
rather than being able to use a react ion- diffusion framework to obtain an understanding of 
the mechanisms that underly the observed macroscopic behavior of Proteus. 

This paper is organized as follows. We provide a brief discussion of Proteus biology, 
followed by a section going into more detail on the critical differences and implications 
between the structured multiscale approach used in this paper and the reaction-diffusion 
models used elsewhere. We state the model equations and their biological meaning. We 
then discuss the numerical methodology used to solve the equations, and present the results 
along with a discussion. We close with our conclusions and suggestions how the validity of 
the underlying assumptions of the structured multiscale models and the reaction-diffusion 
models can be determined empirically. 

2 Proteus mirabilis Swarm- Colony Development 

When inoculated onto an agar surface, Proteus begin formation of strikingly regular spatio- 
temporal patterns that begin with three initial phases: a lag phase, a first swarming phase, 
and a first consolidation phase; these form what appear to be the center circle and first 
ring of a bulls-eye. This is followed by repeating cycles of ring, or terrace, formation that 
consist of a swarming phase followed by a consolidation phase. These subsequent terraces 
are equal in width and time of formation. Perhaps most interestingly, the time of terrace 
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formation is invariant under changes in the glucose or agar concentration of the substrate, 
although changes in agar concentration do affect the ratio of time spent swarming versus 
consolidating, and changes in glucose increase terrace width. Terrace formation does vary 
by temperature [221 ■ 

Varying glucose concentration greatly alters biomass production without strongly affect- 
ing the spatial or temporal aspects of the swarm and the consolidation cycle of terrace 
formation. Changes in agar concentration do not change the total time of terrace forma- 
tion, but higher agar concentrations shorten the time spent swarming and the width of the 
terraces, and lengthen the time spent in consolidation [TU] . 

A short summary of the physiology of Proteus cells can be found in a section similar 
to this one in |2] , and in more detail in 1^ . The critical aspect for understanding the 
macroscopic behavior is the transition between the two cell types: mononuclear cells that 
are motile in a liquid environment and immotile on an agar surface, called "dividing" cells in 
this paper and "swimmer" or "vegetative" cells elsewhere in the literature; and multinuclear 
filament cells which possess longer flagella and are able to move on an agar surface by 
grouping together with other swarmer cells in what are called "rafts" . The process in which 
dividing cells become swarmer cells is called "differentiation." Differentiation only occurs 
above one dividing-cells density and below another. 

When the multinuclear swarmer cells approach a maximum size, they rapidly break down 
into mononuclear dividing cells. This process is called "dedifferentiation." Since mitosis 
results in exponential growth, the size of a swarmer cell is an exponential weight times its 
age. Thus, dedifferentiation occurs when swarmers reach a maximum age. 

We suggest that the differentiation-dedifferentiation cycle of Proteus cells, including a 
density-lag before differentiation occurs and a sharp age near which dedifferentiation occurs, 
is a critical mechanism for generating the macroscopic colony behavior, and thus density- 
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dependent mechanisms alone are insufficient. 

A deeper theoretical and empirical understanding of how deviations from this cell cycle 
affect the macroscopic behavior is needed to resolve the issue. The models presented in this 
paper aim to advance the theoretical understanding, and experiments, perhaps similar to 
those by Hay, Tipper, Gygi, and Hughes piOj where a mutated strain of Proteus was unable 
to swarm due to a lack of sufficiently long flagella, can provide a more sound empirical basis 
for the models. 

3 Critical Differences with Reaction-Diffusion Models 

In this section we discuss the critical difference between the structured multiscale models 
in this paper and react ion- diffusion models of Proteus swarming. We focus on the work by 
Czirok, Matsushita, and Vicsek [Zj. This work makes explicit claims of the sufficiency of 
a reaction-diffusion framework, utilizing only density-dependent mechanisms, to answer all 
questions of interest concerning Proteus swarm-colony development. We direct less attention 
to the paper by Medvedev, Kaper, and Kopell j2lj because their work focused more on the 
interesting mathematical phenomenon of a reaction-diffusion equation with periodic front 
dynamics and did not make such claims about fully understanding the biology of interest. 
The Medvedev, Kaper, and Kopell [21] paper was discussed in more detail in 

The models in [3| consist of a reaction-diffusion equation for the swarmer cells and an 
ordinary differential equation at each point in space for the "vegetative" cells. The term 
"vegetative" cells was initially used in jlj with the same meaning as "dividing" cells in this 
paper, and later used to mean dividing cells and swarmer cells which are unable to contribute 
to swarming. The term "consolidation" has generally been used in the literature to mean the 
macroscopic-scale colony phase when the front is not moving and the newly formed terrace 
is being built up. Czirok, Matsushita, and Vicsek use the term to mean the cellular-scale 
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process of "dedifferentiation" . This use is confusing given the differences in scale, and that 
dedifferentiation involves the breakup of multinuclear filament cells into several mononuclear 
dividing cells. 

The models in [7j use a minimum and maximum population density of dividing cells 
for which differentiation into swarmer cells can occur. The upper limit has been used in all 
models of Proteus swarm-colony development we are aware of to date [21 13 El 1^- Medvedev, 
Kaper, and Kopell first saw the importance of the density lag in differentiation from dividing 
to swarmer cells in the biological literature [211 and incorporated it into a model. The lower 
limit is used here and in [2l[7j, but not in ITU]. One difference between the results presented 
here and in [2] and those in [7j is that Czirok, Matsushita, and Vicsek get a loss of regularity in 
their results if the density window for swarmer cell production becomes too wide. This loss of 
regularity takes the form of the consolidation phase becoming shorter after each cycle. In the 
models here and in the dividing-cell population density is scaled by the density window 
for swarmer cell production in the nondimensionalization. The nondimensionalization as 
described in j2] shows why this is natural. Instead we vary the center of the interval, denoted 
by Vc- As Vc increases, so does total cycle time and consolidation time, while swarm time 
remains invariant. Although we can remove the consolidation phase by changing Vc, the 
result of each simulation retains the temporal and spatial regularity. 

The diffusivity used in [7] is simpler than those used in [TU] or flT^ , but is more complicated 
than the diffusivity used here and in |2]. 

The models in this paper extend those developed in to examine the necessity of a 
sharp age of dedifferentiation from swarmer to dividing cells. The critical difference between 
the models in this paper and in |2| and those in [7j is the mechanism of dedifferentiation. 
We argue that the assumptions underlying the models are mutually exclusive. The models 
in this paper require that dedifferentiation occur almost exclusively near a critical age. The 
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models in [7j use a constant rate of dedifferentiation. 

The important differences and commonalities between the age explicit frameworks used 
here and those of Esipov and Shapiro jlO|| were discussed in much detail in A crucial 
difference relevant to this discussion is that we dispense with a memory field in the diffusion, 
which constitutes a built in hysteresis. Unfortunately, Czirok, Matsushita, and Vicsek [7] 
focus on the memory field, rather than the critical insight of Esipov and Shapiro of including 
the cell cycle in the model, in their discussion of ,10.|. The major deficiencies of the models 
and computations in PU] are covered in much detail in However, in the context of [7j, 
this paper and [2^ constitute a defense of the need to explicitly represent age and the cell 
cycle in a model of Proteus swarm-colony development. 

The models in j2] reproduce three major aspects of Proteus swarm-colony development: 
temporal regularity of overall cycle time, spatial regularity, and control of the ratio of swarm 
time to consolidation time. By retaining the explicit age structure in the mathematics, the 
models reproduce the colony behavior by relying on two main mechanisms: a density-lag in 
differentiation from dividing cells to swarmer cells, and a sharp age of dedifferentiation from 
swarmer cells to dividing cells. The react ion- diffusion models of both in Czirok, Matsushita, 
and Vicsek jTj and Medvedev, Kaper, and Kopell |2Ij reproduce only the temporal regularity 
of the overall cycle time and the spatial regularity. No attempt is made to address the 
ratio of swarm time to consolidation time. Given that this is controlled in the models 
in j2] by an age threshold, it is not clear how such a control can be inserted into one of 
the reaction-diffusion modeling frameworks. However, we concede that a control within a 
react ion- diffusion framework is within the realm of the possible. 

Instead we argue the necessity of including the cell cycle explicitly by showing that 
moving from a sharp age of dedifferentiation to a broader distribution of dedifferentiation 
ages causes a breakdown in regularity, which would exclude density-dependent mechanisms 
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alone from explaining Proteus swarm-colony and indicate the need to include age explicitly. 
The models would be verified if experiments with strains of Proteus which do not have a 
sharp age of dedifferentiation give rise to irregular colony formation. 

4 The Models 

The structured multiscale models in this paper are based on a mathematical framework 
with some history behind it. At the colony or population scale, Skellam |21] considered the 
effects of diffusion in his classic work of 1951. At the individual scale, Sharpe and Lotka in 
1911 and McKendrick in 1926 considered population models with linear age structure |25j . 
More recently, Gurtin and MacCamy jT2I considered models with nonlinear age structure. 
Rotenberg j2S] and Gurtin [TT] posed models dependent on both age and space. Gurtin and 
MacCamy 13^ differentiated between two kinds of diffusion in these models: diffusion due to 
random dispersal, and diffusion toward an area of less crowding. Existence and uniqueness 
results can be found for various forms of these models in Busenberg and lannelli [SJ, di Blasio 
jSj, di Blasio and Lamberti Langlais and MacCamy ^H]- Further analysis has been 
done by several authors [13 CEl UHl 1201 • 

We present in this section dimensionless equations, variables, and parameters. We assume 
the colony is radially symmetric and scale space so that radius goes from zero to one. We 
scale time and age by the time it takes a cell to subdivide (typically 1.5 hours.) A complete 
discussion of the nondimensionalization is contained in [2]. The dimensionless variables r, a, 
and t represent radius in two-dimensional space, age, and time, respectively. The function 
u{r,a,t) represents the swarmer-cell population density at radius r, age a, and time t. The 
functions p{r,t) and v{r,t) represent the biomass density of swarmer cells and the dividing- 
cell population density, respectively, at radius r and time t. 
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The nondimensional model consists of the age-structured nonhnear diffusion system 

dtu + daU = ^dr{rdr{D{p)u)) - n{a)u, < r < 1, a > 0, t > 0, (4.1a) 

POO 

dtV = {l-i{v))v+ I n{a)ue'' da,0 <r <l,t> 0, (4.1b) 

Jo 

u{r,0,t) = C{v)v{r,t),0<r <1, t>0 (4.1c) 

with boundary and initial conditions 

dr{D{p{l,t))u{l,a,t)) ^0, a>0,t>0, (4.1d) 
u{r, a, 0) = 0, < r < 1, a > 0, (4.1e) 
v{r,0)^vo{r), < r < 1. (4.1f) 

The total motile swarmer cell biomass is given by 

u{r, a, i)e" da, <r <1, t>0. (4.1g) 

min 

Proteus move through a process of raft building that requires two things: sufficient 
maturity in swarmer cells to contribute to raft building, and a sufficient biomass of mature 
cells to form the rafts. The lower limit of integration, a^i„, is the minimum age when a 
swarmer cell is sufficiently large, with sufficiently long flagella, to contribute to motion on 
the agar surface. This parameter controls the ratio of swarm time to consolidation time 
without changing the total time of ring formation. The parameter a^^^ is related to agar 
concentration. The higher the concentration, the drier the surface, raising the value of a^^^. 
Higher agar concentration, and thus higher a^i„, shortens the swarming phase and lengthens 
the consolidation phase in terrace formation. 

The diffusivity has the form 

D{r,t) = L'omax{(p(r,i) -p^i„),0}. (4.1h) 



The parameter p^^^ is the minimum biomass needed for swarmers to build rafts capable of 
moving on the agar surface. 

We use a differentiation function with a lag phase that is a piecewise cubic with 
support of length 2: 



The interval [fc — l,fc + 1] is the swarmer-cell production window. The interval width of 
2 was obtained in the nondimensionalization. The parameter Vc represents the lag phase 
in swarmer-cell production [21]. We use a compactly supported cubic function because it 
is smooth. The important aspects of the functional form of ^ are that it is zero above and 
below certain thresholds. An appendix in ^ shows that the nature of the results do not 
depend strongly on the shape of the curve. The shape of the curve does help with numerical 
issues concerning the degenerate diffusion. The parameter has no analog in |10j . 
We take the initial condition to be 



Again, it is mass, not shape, that matters qualitatively. Shape effects numerical efficiency. 

The major difference in the models in this paper and those in j2] is the dedifferentiation 
modulus, and the resulting behavior of the system to changes in dedifferentiation from 
swarmer to dividing cells. In 0, we used a form that was the sum of a Heaviside function 
and a Gaussian distribution so that in the limit of the spread parameter cr ^ it would be 
transparent that we obtain a model with no explicit function /i but could instead represent 
the sharp age of dedifferentiation by an upper limit in the age domain. Only the limiting 
case was treated in (2jj. 

Our goal in this paper is to show the necessity of retaining explicit space structure and 
a sharp age of dedifferentiation. Consequently, we examine the change in swarm-colony 




|f — fc| < 1 

otherwise. 
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regularity as dedifferentiation become less sharp. We use a function with similar shape but 
a simpler expression than the dedifferentiation modulus used in P] , 



where a is the spread parameter and /io is a height parameter. This function is equivalent 
to a Fermi distribution. For a fixed /io and taking the limit cr — ^ 0, we get a situation 
where swarmer cells dedifferentiate into their component dividing cells all at the same age, 
a^iax- For a small but not equal to zero, this dedifferentiation modulus represents a situation 
where the probability of dedifferentiation is low for young swarmers, increases rapidly as 
they approach a^^x in age, and remains high afterwords. What is important is the rapid 
change in the dedifferentiation rate. Simulations carried out with a piecewise continuous 
linear function /i, that is zero below a^^^ — a, /xq above a^^x + o^, and increases linearly in 
between, show no qualitative and minor quantitative differences with the form in equation 



It may be unclear at first why the magnitude coefficient of /x (the ratio fio/cr) should 
increase as the hyperbolic tangent approaches a step function. The intent is to model the 
effect on colony behavior of moving away or toward a sharp age of dedifferentiation in the 
model. We show that if we keep /io/cr fixed as we vary a, we converge to something other 
than a sharp age of dedifferentiation as a — >■ 0. The probability that a swarmer cell has not 
dedifferentiated by age a is given by 




(4.1k) 





(see pg. 82 of [H].) Using the specific form of fi in equation ()4.1k|) . we obtain 




(4.3) 
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Figure 1: The probability n(a) that a swarmer cell has not dedifferentiated by age a. In 
figure [T(a)j /io = 1-0 in equation ()4.H|1 . In figure [T(bj| /iq = a in equation ()4.Hj) . These figures 
illustrate the need to increase the magnitude of the dedifferentiation modulus to approximate 
a sharp age of dedifferentiation from swarmer cell to dividing cell in the models. 

Figure H illustrates how keeping /xq fixed as a ^ yields convergence of 11 to a step function, 
whereas keeping /io/o" fixed gives convergence to a function with decay to zero rather than 
a discontinuity. 

The correct way to represent a sharp age of dedifferentiation has proven to be nontrivial. 
Models in (specifically "Model A" as defined on pp. 252-253) conflate the desired 
probability density function (pdf) for dedifferentiation ages (a Dirac delta function centered 
at the sharp age of dedifferentiation) with the dedifferentiation modulus. As can be seen by 
setting ^{oj) = 6{uj — a,-,-,^^) in equation 14. 2[ a Dirac delta function as the dedifferentiation 
modulus would give rise to a situation where there is no dedifferentiation until a^^x, at which 
point the population is reduced by a factor of 1/e. This remaining swarmer-cell population 



12 



never dedifferentiates, and continues to age (and hence grow), resulting in no consolidation 
period whatsoever. 

Esipov and Shapiro PU] attempted to examine the importance of a sharp age of dedif- 
ferentiation by treating two models, "Model A" with dedifferentiation as mentioned above, 
and a "Model B" with a constant dedifferentiation modulus. As discussed in P], "Model 
A" , when solved accurately, does not reproduce the observed Proteus colony behavior, which 
limits the utility of that treatment. Just as important, only considering the extreme cases 
leaves open several questions. All swarmer cells do not dedifferentiate at precisely the same 
age; this is a construct useful for modeling. Does even a minor deviation from a sharp age of 
dedifferentiation lead to colony irregularity? If so, this would invalidate the model. Also left 
open is the nature of the transition from regularity to irregularity of colony formation. We 
attempt to answer these question using mechanisms in our dedifferentiation function absent 
from those in |lUj . 

5 Computation and Results 

In this section we briefly describe the computational methods used to obtain our solutions 
and then present those solutions with their biological interpretations. 

The computational methods and software used in this paper are essentially those used to 
solve the model equations in The main exception is the use of a modulus, /i, rather than 
a finite age domain truncated at the septation age, to model dedifferentiation from swarmer 
to dividing cells. This required changes in how information was passed from the dividing- 
cell equation to the swarmer-cell equation. Previously, the population of cells undergoing 
dedifferentiation could be determined from the oldest age cohort at each point in space at 
the end of the time step and passed to the dividing-cell equation at the start of the next 
time step. For the models in this paper, the population of dedifferentiating swarmer cells 

13 



had to be computed across all age cohorts at each point in space and then aggregated for 
passing to the dividing-cell equation. 

A more detailed description of the numerical methodology as specifically applied to Pro- 
teus can be found in 0, and the original formulations and analyses can be found in P1|3J|3]. 
To summarize, the software uses a moving-grid Galerkin method in the age variable to de- 
couple age from time, which results in a degenerate nonlinear parabolic partial differential 
equations in space and time for each age cohort. The parabolic system is then solved using 
a step-doubling extrapolation method jl]. 

We emphasize the importance of reliable numerics. In the systems presented in this paper 
and in 0, the regularity in terrace formation is due to the regularity in the aging of the 
swarmer cells. A coarse, stage-structured numerical approximation to the continuous aging 
term, such as that used in [lOj, can constitute a qualitatively different swarmer-cell life cycle 
than that of the original continuous model. Also, the degenerate diffusion can be altered 
by the choice of the time step. The numerical computation may then show periodic front 
dynamics that are not obtained by an accurate solution of the original continuous equations, 
but are induced by a regularity in the numerics. An appendix in |2] details how accurate 
computations change the behavior of the Esipov-Shapiro system from what is presented in 

uni. 

A uniform spatial discretization of size Ax = 1/300, and a uniform age discretization of 
Aa = 1/40, using piecewise constant basis functions, was sufficient to solve the system within 
a relative error in the L^-norm of less than 1%. A uniform age discretization in the context 
of a moving grid method means that all but the first and last age intervals are constant in 
length, and that a new age interval is introduced at the birth boundary when the old birth 
interval reaches Aa in length. Convergence in time was obtained by adjusting a tolerance 
parameter for the adaptive time-stepping in the step-doubling algorithm so that the relative 
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error in the L^-norm was also less than 1%. 

The goal of this paper is to model the effects of a non-sharp dedifferentiation age on colony 
regularity. Consequently, we focus on the behavior of the system defined by equations ()4.1ap - 
fl4.1k|) as we vary fiQ and a. For the computations in this paper, we set a„,^x = 2.67, Vc = 8.0, 
Do = 2 X 10~^, p^in = 0.5, ^0 = 0.5, and = 0. 

The response of the system to changes in these six parameters was investigated in detail 
in j2] (we have dropped the "hat" or capital notation in this paper, which was used in 
[2] to denote the nondimensional version of a variable or parameter.) In summary, a^^^ 
is the nondimensional sharp age of dedifferentiation (in this paper its meaning changes to 
become the "center" of a mollified step function.) As a,^^^ increases, so does the time spent 
swarming, the terrace width, and the overall cycle time. The time spent in consolidation 
decreases. The parameter Vc is the nondimensionalized lag in swarmer cell production, ^. 
As Vc increases, consolidation time, total cycle time, and terrace width increase, where as 
swarming time remains nearly invariant. The constant of diffusion, Dq, represents changes 
in glucose concentration and has a significant effect only on terrace width. The parameter 
is the constant of the differentiation ratio ^, and lies between zero and one. As increases, 
so does swarming time and terrace width, whereas total cycle time and consolidation time 
decrease. 

The parameter a^j^ is the minimum age at which a swarmer cell can contribute to raft 
building and movement and controls the ratio of swarm time to consolidation time within a 
cycle. It only has meaning in the context of an age-explicit model. As agar concentration 
increases, so does the difficulty of moving on the substrate, which is represented by increas- 
ing a^i^. As a^in increases, swarming time and terrace width decrease, consolidation time 
increases, and total cycle time remains the same. 

The parameter p^^^ is the minimum swarmer biomass needed for raft formation and the 
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onset of swarming. It is thus measured in the same units as population density, and differs in 
an important respect from a,„in, which is measure in the same units as age. This difference is 
worth noting given that we argue the necessity of including both density- dependent processes 
and the cell cycle in a model, rather than relying on density-dependent mechanisms alone to 
explain the observed Proteus behavior. The response of the system to changes in p^^^ differs 
from changes to a^^^. As p^^^ increases, swarming time, total cycle time, and terrace width 
decrease, whereas consolidation time increases. The parameter p^^^ is needed to create a 
distinct consohdation phase. 

To understand the breakdown of colony regularity as the age of dedifferentiation become 
less sharp, we examine first the response to changes in ^q. Decreasing for a fixed a 
has the effect of decreasing the amount of dedifferentiation after a^^^, which can ultimately 
result in a continuation of swarming and a loss of a consolidation phase. As shown in 
figure this breakdown can be quite irregular. The temporal breakdown occurs before an 
apparent breakdown in spatial regularity. Not shown in the figure is the effect of increasing 
Hq. Increasing /io above what is shown results first in regular colonies with narrow terrace 
widths, and then the loss of regularity through the loss of any meaningful swarming phase. 
This is intuitive since increasing /io results in earlier and more frequent dedifferentiation, 
which in turn results in a lack of large swarmer cells to contribute to raft building and the 
consequent swarming. 

We next examine the response to changes in a. This parameter study corresponds to 
the change in the probability function n(a) that a swarmer cell has not dedifferentiated by 
age a shown in figure |l(a)} As illustrated in figure El colony regularity breaks down as the 
age of dedifferentiation becomes less sharp. It is important to note that the colony remains 
regular within an interval around a = 0. As mentioned above, the swarming cycle becoming 
irregular for a near zero, but not equal to zero, would invalidate the model. It is interesting 



16 



(a) Simulated Colony 




Time 



(b) Position of Colony Front 



Figure 2: Computed Proteus mirahilis swarm colonies showing a breakdown of spatial and 
temporal regularity as the dedifferentiation function becomes less sharp. Note that temporal 
regularity breaks down before the apparent spatial regularity. Figure 2(a) consists of the 3D 
plots of the swarmer- and dividing-cell biomass viewed from directly above with a combina- 
tion of ambient light, diffuse reflection, and specular reflection. The time of the snapshot 
is the earlier of t = 30 or when it hits the boundary. Figure |2(b)| shows the corresponding 
radius of the swarm colonies as functions of time. Areas where the function has zero slope 
correspond to consolidation periods in the swarm-colony development. The dedifferentia- 
tion parameter yUo varies for the three colonies, left to right, from /io = 0.1 to yUo = 0.01 
to /io = 0.001. Parameters common to all three swarm colonies are a = 0.01, a^^^^ = 2.67, 
Vc = 8.0, Do = 2 X 10~^, p^i,, = 0.5, = 0.5, and a^i^ = 0. 
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Figure 3: Swarm-colony front dynamics showing the breakdown of spatial and temporal 
regularity as the dedifferentiation age becomes less sharp. Parameters common to all three 
swarm colonies are /i = 1.0, a^^^ = 2.67, Vc = 8.0, Dq = 2 x 10~^, p^i^ = 0.5, = 0.5, and 

C^min 0. 

to note that the breakdown of colony regularity is nearly hysteretic as a function of a. In the 
results shown in figure El we set /iq = 1. Results do not change quahtatively for /io near this 
choice, however, as shown in figure |2l /io is an important parameter in determining colony 
regularity. Fortunately, we have the shape of 11 (a) as a means to determine which parameter 
choices represent a sharp versus a mollified age of dedifferentiation. 

The last parameter exploration we present is that of fixing the ratio fj^o/cr and varying 
a. This corresponds to the change in the probability function 11(a) shown in figure [T(b)] As 
0" — ^ 0, the colony formation never becomes truly regular. This irregularity corresponds to 
a lack of sharpness in n(a) as cr ^ 0. Figure IH illustrates our results for /io = cr. Setting 
/io = 3cr results in a loss of a swarming phase as a increases rather than the loss of a 
consolidation phase. 
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Figure 4: Swarm-colony front dynamics showing irregularity of all colonies when /io/cr is 
fixed. For the colonies shown in this figure, /iq = cr. Parameters common to all three swarm 
colonies are a^^^ = 2.67, Vc = 8.0, Dq = 2 x 10~^, p^^^ = 0.5, = 0.5, and a^i^ = 0. 

6 Conclusions 

This paper constitutes an argument that the regularity in the Proteus cell cycle of differ- 
entiation and dedifferentiation between immotile dividing cells and motile filament swarmer 
cells, along with mechanisms that depend on cell density, underlies the temporal and spatial 
regularity in swarm-colony development. This is at odds with the implications of reaction- 
diffusion models which attribute Proteus behavior to density-dependent mechanisms alone. 
The first conclusion to make is that the question of the correct modeling approach might be 
answered empirically by studying the colony regularity of strains of Proteus that do not have 
a sharp age of dedifferentiation. The probability that a swarmer cell has not dedifferentiated 
by age a, denoted by n(a), can be used as a guide to determine experimentally the forms of 
the dedifferentiation modulus /i in equation ()4.1ap . 

Density-dependent effects do matter in understanding Proteus swarm-colony develop- 
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ment. Experiments to determine how differentiation from dividing cells to swarmer cells 
depends on local population density would be particularly valuable in finding the forms of 
the differentiation function, ^. 

We, along with the other scientists modeling Proteus swarm-colony development of whom 
we are aware, used diffusion to model movement. Diffusion may not be the best represen- 
tation of the spatial dynamics. Areas of future research are the specifics of how the much 
larger and broader issue of biological motion relates to Proteus motion on an agar surface, 
and extending the models to explicit two-dimensional space to understand the interaction of 
multiple colonies. 

We close by noting that the likely importance of the cell cycle in determining the macro- 
scopic behavior makes Proteus a model system for a multiscale understanding of other bio- 
logical systems where the spatial or temporal behavior is a manifestation of local kinetics. 
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